Radiative Scalar Meson Decays in the Light-Front Quark Model 
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We construct a relativistic 3 Po wavefunction for scalar mesons within the framework of light- 
front quark model(LFQM). This scalar wavefunction is used to perform relativistic calculations of 
absolute widths for the radiative decay processes (0 ++ ) — » 77, (0 ++ ) — > 07, and (0 ++ ) — » p-y 
which incorporate the effects of glueball-gg mixing. The mixed physical states are assumed to 
be /o(1370), /o(1500), and /o(1710) for which the flavor-glue content is taken from the mixing 
calculations of other works. Since experimental data for these processes are poor, our results are 
compared with those of a recent non-relativistic model calculation. We find that while the relativistic 
corrections introduced by the LFQM reduce the magnitudes of the decay widths by 50-70%, the 
relative strengths between different decay processes are fairly well preserved. We also calculate 
decay widths for the processes cj> — ■» (0 ++ )7 and (0 ++ ) — > 77 involving the light scalars /o(980) and 
ao(980) to test the simple qq model of these mesons. Our results of qq model for these processes 
are not quite consistent with well-established data, further supporting the idea that /o(980) and 
ao(980) are not conventional qq states. 



I. INTRODUCTION 

As is well known, the assignment of the scalar ( J PC = 
++ ) qq states has long been an enigma in hadron 
spectroscopy. Unlike the elegant, ideally mixed vec- 
tor and tensor multiplets, it is still controversial which 
are the members of the expected L = S = 1 qq mul- 
tiplet since there are now too many ++ mesons ob- 
served in the region below 2 GeV for them all to be 
explained naturally within a qq picture^. For ex- 
ample, 2 isovector(/J PC ' = 1 ++ )[a (980), a (1450)] 
and 5 isoscalar(0 0++)[/ (600)(or a), /o(980), /o(1370), 
/o(1500), /o(1710)] states have been reported by the Par- 
ticle Data Group lj. This has led to the suggestion that 
not all of them are qq states. The main reason for this 
situation is that around the relevant mass region there ex- 
ist other kinds of particles such as KK molecules 0,0,0; 
glueballspllfl, four- quark (qqqq) systems)?}, or hybrids. 

Interpreting the structure of each of the known scalars 
has proven to be a fairly controversial endeavor. Take, 
for example, the light scalars (i.e. those below 1 GeV). 
Due to some of the difficulties associated with /q(980) 
and do (980) — e.g. the strong couplings to KK in spite 
of their masses being at the KK threshold, and the 
large discrepancies of 7T7t0|,770 and </>-radiative decayQ 
widths between non-relativistic (NR) quark model pre- 
dictions and experimental data — Weinstein and IsgurQ 
proposed the isospin 7 = /o(980) and the 7=1 ao(980) 
within the NR potential model as the "KK molecules" . 
However, a more conventional interpretation for these 
states has been given by Tornqvist and Roos0 who 
analyzed the data on the /o(980), / (1370), a (980) and 
ao(1450) as unitarized remnants of qq 1 3 Pq states with 
six parameters and theoretical constraints including fla- 
vor symmetry, the OZI rule, the equal-spacing rule for the 
bare qq states, unitarity, and analyticity. In this work, 
the authors concluded that the /o(980) and the /o(1370) 
are two manifestations of the same ss, while the ao (980) 



and the eto(1450) are two manifcstaions of the same ud 
state. 

The interpretation of the structures of some of the 
heavier scalars has been somewhat less controversial. For 
example, /o(1370) and ao(1450) are most commonly in- 
terpreted as qq states, even though the flavor assignments 
for each are still unclear. Also, there appears to be gen- 
eral agreement on identifying /q(1500) as a glueball, pos- 
sibly mixed with nn = (uu + dd) j y/2 and ss. This inter- 
pretation has followed both from lattice QCD which, in 
the quenched approximation, predicts that the lightest 
glueball has J PC = 0++ and a mass of 1.55-1.74 GeV 
[lfjl PL1|| . as well as from the fact that /o(1500) decays 
strongly into 7T7T but not into KK. 

Most recently, Close and Tornqvist ^1 have proposed 
a scheme which sorts the light scalars into two distinct 
nonets: one nonet above 1 GeV and another below 1 
GeV, with different physics operating in each. The nonet 
above the 1 GeV threshold is comprised of the qq states 
mixed with the scalar glueball. The glueball's presence is 
inferred from the overpopulation of isoscalars in this mass 
region. The nonet below 1 GeV is made up of qqqq and 
meson-meson molecules. As such, /o(980) and ao(980) 
can be thought of as superpositions of four-quark states 
and KK molecules. The authors of Ref. ^3 demonstrate 
that such a scheme involving two scalar nonets can be 
described using two coupled linear sigma models. 

In this paper, we investigate various radiative scalar 
meson decays which could provide important clues on the 
internal structures of these states |a.ll3t- For the calcula- 
tions, we employ the light-front quark model (LFQM) 
which has been used very successfully in the past to 
compute decay rates for pseudoscalar, vector, and axial- 
vector mesons 0, 0, ^3 ■ Extending the model to in- 
clude scalar states mainly involves the construction of a 
new 3 Po light-front model wavefunction. In section II, 
we give the detailed form of this wavefunction and ex- 
plain how the model parameters are obtained. This new 
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scalar wavefvmction is used to study the radiative decays 
of the heavy scalars / (1370), / (1500), and / (1710), as 
well as the light scalars ao(980) and /o(980). 

In the case of the heavy scalars, we adopt the scheme 
of Close and Tornqvist in which /o(1370), /o(1500), and 
Jo (1710) are considered to be mixtures of nn, ss, and 
gg. The flavor-glue content of each state is taken from 
mixing analyses done by Lee and Weingarten |17| and by 
Close and Kirk ^^|. Taken together, these works provide 
mixing amplitudes for three distinct cases of the scalar 
glueball mass: (1) a heavy glueball (M nn < M sS < M gg ), 

(2) a medium weight glueball {M nfl < M gg < M sS ), and 

(3) a light glueball (M gg < M nfL < M sS ). The details of 
these three mixing scenarios are outlined in section III. In 
the case of the light scalars, ao(980) and /o(980) are as- 
sumed to be conventional qq states. The flavor content of 
ao (980) is then (uu — dd) / y2, and that of /o (980) is some 
superposition of nn and ss. Rather than attempting to 
determine the degree of mixing for / (980), it suffices to 
examine the two ideally mixed cases: Jo (980) = nn and 
/ (980) = S s. 

In section IV, the general forms of the Q 2 -dependent 
transition form factors for the processes (1 ) — > 
(0++) 7 *, (0++) -> (I -- h*, and (0++) -> 77* are de- 
rived. In the limit as Q 2 — > 0, these form factors yield 
the decay constants for the real photon processes which 
can then be used to compute the corresponding decay 
widths. In section V-A, we present our numerical results 
for the heavy scalars. This includes the form factors and 
decay widths for the specific processes f — > (fij, / — > pj, 
and /o — > 77. Our results for the light scalars involved 
in the processes 4> — * /o(<Zo)7 and fo( a o) —> 77 are given 
in section V-B. A summary of the paper's salient points 
and a brief discussion of our intended future work is given 
in section VI. In the Appendix, the explicit form of the 
trace used in section IV. A is presented. 



II. THE MODEL WAVEFUNCTIONS 

One of the popular quark models in the light-front for- 
malism is the invariant meson mass(IM) scheme^l |2(j 
in which the invariant meson mass square Mq is given by 



2 k 2 

*V* I 



E 



(i) 



In our analysis, we will only consider the light- meson 
sector(w, d, and s quarks) with equal quark and anti- 
quark masses (m q — m q ). 

The light-front qq bound-state wavefunction of the 
scalar( 3 Po) and vector( 3 Si) mesons can be written in the 
following covariant form 

*M(^,k,-i, A*) = u\ q {p q )T M vx i {pq)(j)M{x l ,'k il _) 

= '^-\ q \ 9 ( x i,' k i±) ( l ) M(xi,'k l± ), (2) 

where IZf 1 \_ is the spin-orbit wavefunction, which is ob- 
tained by the interaction independent Melosh transfor- 



mation from the ordinary equal-time static spin-orbit 
wavefunction, and 4>M(xi,'k.i±) is the radial wavefunc- 
tion. The operators Tm for the scalar(S) and vector(V) 
mesons are given by 



r s = 



(M$ + m q + m ? )^/2[(M s )2 - (m 2 - mf)] 

z (Fv+M[) Ah) 

{M% + mq + m g )^2[{MX) 2 -(ml~m\)] 



,(3) 



where P(s,v) = (Pq +Pq), K = {Pq — Pq)/^ is the relative 
four-momentum between the quark and antiquark, and e 
is the polarization four-vector of the vector meson(with 
momentum Py), which is given by 



e"(±) 
ex(±) 
c"(0) 



= T 



e , e ,ej_J 
(l,±i) 



0,- T e ± (±)-P v ±,e±(±) 



V2 



1 

My 



Ml 



V ! ' 



1 T 



(4) 



The operator Ty was derived in Ref. 1211 . We followed 
the same procedure, detailed in Ref. [2l|, in order to 
derive Tg. Note that in the case of the vector me- 
son, the operator T v has the expected form, (P + M) tf. 
However, because the scalar ( 3 Po) state possesses non- 
zero orbital angular momentum, the proper Melosh- 
transformed spin-orbit wavefunction is not simply given 
by Ts = (P + M), as one might expect. The form is more 
complicated as shown in Eq.Q, and it depends explic- 
itly upon the relative momentum between the meson's 
constituents. This same type of behavior was demon- 
strated for the axial-vector meson in Ref. [21| ■ Since the 
axial-vector ( 3 Pi) state also possesses non-zero orbital 
angular momentum, the spin-orbit wavefunction is not 
simply given by (P + M) ^y 5 . The correct form contains 
an additional factor which explicitly depends on the rel- 
ative momentum between the quark and anti-quark. It 
is interesting to note, however, that in the case where 
m q = m q = m (which we use throughout this work), the 
expressions in Eq. © reduce to 



2V2A/ S 



T v = 



-1 



\ + MX)Ah). (5) 



So, in the equal mass case, T s does have the expected 
form. We confirmed the similar reduction of the axial- 
vector meson wavefunction in the equal mass case. These 
expressions can be further simplified to the form we will 
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use in our analysis: 
1 



r s = 
rv = 



2^2' 
-1 



e ' (Pg - Pg) 



2m 



(6) 



The spin-orbit wavefunctions satisfy the following rela- 
tions 



^ -V ^^Ag Ag^'Ag Ag 
Aq Aq 



l(M 2 -4m 2 ) = |k| 2 , 



/ Ag Ag 



Aq Aq 



= 1, 



(7) 



where k = (kx,/c z ) is the three momentum of the con- 
stituent quark and k z — (x — \)Mq. Note that the total 
wavefunction ^ g(x, k^) for the scalar meson vanishes at 
|k| = in accordance with the property of P-wave func- 
tion. 

For the radial wavefunctions 4>m(xi, kij_), we shall use 
the following gaussian wavefunctions for the scalar and 
vector mesons 



(f>s(x,'k A 



dk z 



= ^^exp(-k 2 /2/3 2 ), (8) 

where J\f = 4(^-) 3 / 4 and dk z /dx is the Jacobian of the 
variable transformation {a;,kx} — > k = (fc z ,kj_) defined 
by 



9^ 
dx 



A/n 



4x(l - x) ' 



(9) 



and the normalization factors are obtained from the fol- 
lowing normalization of the total wavefunction, 



fix 



d 2 k ± 
16tt 3 



|* M (^,k ± )| 2 = 1. 



(10) 



The wavefunctions depend on only two model param- 
eters: the constituent quark mass, m, and the binding 
strength, (3. For these parameters, we use the values de- 
termined in Ref. 0] which are (m Uj d = m n = 0.22 GeV, 
u4 = f3 n = 0.3659 GeV) and (m s = 0.45 GeV,/3 s = 
0.4128 GeV). These values were obtained by fitting the 
LFQM spectrum for pseudoscalar and vector mesons — 
obtained using a QCD-inspired model Hamiltonian with 
a linear confining potential — to experimental data. In 
this paper, we will use these values of the parameters for 
both vector and scalar mesons. In a future work, we in- 
tend to perform a separate fit to scalar meson data to 
see whether or not the model parameters would differ 
significantly. 

Since our analysis deals with <f> decays, a value for the 
uj-<j> mixing angle is necessary. We use 8^-4, — ±7.8° 



(the sign cannot be fixed), which was also determined in 
Ref. [l4(. This value was obtained by a mass squared 
mixing analysis in which it was assumed that 

I <t>) = - sin r5 w _0 1 nn) - cos S^.^ \ ss) 
\cj) = cos Su-^nn) - sin S^.^ss) , 

and in which the masses of the bare quarkonia were de- 
termined using the model parameters and Hamiltonian 
mentioned in the previous paragraph. 



III. SCALAR MIXING AMPLITUDES 

Glueball-gg mixing can be described using a mass mix- 
ing matrix. Written in the \gg), \ss), \nn) basis, this takes 
the form [H El 



M = 



M gg f V2fr 

f M ss 
V2fr M nn 



(11) 



where / = (gg\M\ss) and r = (gg\M\nn) / y2{gg\M\ss) . 
As described in the introduction, the mixing is assumed 
to be among the three isoscalar states /o(1370), /o(1500) 
and / (1710). These mixed physical states can be written 
as 



|/ (1710)> \ 
l/o(1500)> 
l/o(1370)> / 



a\ b\ c\ 

02 &2 C2 

03 b 3 C 3 




\99) 

ss) I , (12) 



The mixing amplitudes, a*, bi, and Cj can be written 
in terms of the physical masses (M/.), the bare masses 
(Mgg, M sS , and M nfi ), and the glueball-bare quarkonia 
mixing strengths (/ and r). For the present analysis we 
will adopt the values obtained by Lee and Weingarten, 
and by Close and Kirk. 

Beginning with the lattice QCD-motivated assump- 
tion that the bare scalar ss is lighter than the scalar 
glueball (i.e. M nn < M sS < M gg ), Lee and Weingarten 
obtained the following mixing amplitudes |17| : 



0.86 ±0.05 
-0.13 ±0.05 
-0.50 ±0.12 



0.30 ±0.05 
0.91 ±0.04 
0.29 ±0.09 



0.41 ±0.09 
-0.40 ± 0.11 
0.82 ± 0.09 



(13) 



Using a different approach, Close and Kirk jT^| examined 
the constraints placed on the flavor content of /o(1370), 
/o(1500), and /o(1710) by decay branching ratios to pairs 
of pseudoscalar mesons. From a % 2 analysis of the avail- 
able branching ratio data, they obtained various solutions 
depending on which parameters were left free at the out- 
set. One solution was consistent with a glueball which 



lies just above the bare nn (i.e 
It's associated mixing amplitudes are 



M nn < 



Mgg < M S s). 



0.39 ±0.03 
-0.65 ±0.04 
-0.69 ±0.07 



0.91 ±0.02 
0.33 ±0.04 
0.15 ±0.01 



0.15 ±0.02 
-0.70 ±0.07 
0.70 ±0.07 



(14) 
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Another of their solutions was consistent with an even 
lighter glueball which lies below the mass of the nn (i.e. 
M gg < M nn < M sS ). The mixing amplitudes for this 
solution are 



0.25 
-0.37 
-0.89 



0.96 0.10 
0.13 -0.92 
0.14 0.44 



(15) 



In the next section, we shall obtain expressions for 
the transition form factors for various radiative decays. 
Then, in the following section, we will use the mixing 
amplitudes given in Eqs. (|13f) -H15 [) above to numerically 
evaluate the form factors and corresponding decay widths 
for the heavy, medium, and light weight glueball cases re- 
spectively. 



q 



Pi 




Po 



FIG. 1: Primary diagram for A(Pi) -> X(P 2 ) +7(9). There 
is an additional diagram in which the virtual photon interacts 
with the antiquark. 



Our analysis is based on the standard light-front 
frame(g+ =0)13. 



IV. FORM FACTORS FOR RADIATIVE 
DECAYS 

A. The process V(S) ~* S(V) + 7 

The coupling constant g A x-y for the radiative 
Pa (gig) — > Px(g2g)7 decays between vector (V) and 
scalar (S) mesons, i.e. [A,X) = ( 3 Si, 3 P ) or ( 3 P , 3 Si), 
is obtained by the matrix element of the electromagnetic 
current J M which is defined by 

Mi = (X(P 2 )|e 7 - J\A(Pt)) 

= eg A x~{[(e-f ■ ey)(Pi • q) - (e 7 • Pi)(ey • g)], 

(16) 

where e 7 and ey are the polarization vectors of the pho- 
ton and the vector meson, respectively. Since the J z = 
state of the vector meson cannot convert into a real pho- 
ton, the ey should be transversely polarized( J z = ±1) to 
extract the coupling constant <?ax 7 - In other words, the 
possible helicity combinations in the transition V(S) — > 



S(V)"f are either from (J2 



+1,J? 



1 



(J2 — l) JY 

is given by H 



-1). The decay width for A 



T(A^X + j) = a 



_9axj_ 
2J A + 1 



Ml 



2M A 



or from 



(17) 



where a is the fine structure constant and J a is the total 
angular momentum of the initial particle. 

In LFQM calculations, we shall analyze the virtual 
photon(7*) decay process so that we calculate the mo- 
mentum dependent transition form factor, F A x-y*(Q 2 )- 
The coupling constant, g A x~j, can then be determined 
in the limit as Q 2 — > (i.e., g A xy — Pax 7 *(Q 2 = 0)). 
Figure ^ shows the primary Feynman diagram for this 
process. The amplitude is 



Mt = (X(P 2 )\J^\A(P 1 )) 
= ePAx 7 (g 2 )[4(Pi-g) 



P?(evq)}. (18) 



Pi 

Q 
Pi 



Ml 



(P 1 +Pf,P 1 x) = (Pi + ,-^,O x ) 



P, H 



(0, 
Pi 



Ml 



Ml 



Q 2 



q 



p+ 



Q 2 



Pf 



-qxj 



(19) 



-g 2 is the space- like photon momen- 



where = Q 2 
turn transfer. 

The quark momentum variables in the q + = frame 
are given by 



pi 


= (1 


-x)P+, 


p+ = xP+, 


pi± 


= (1 


- x)Pi± 




pi 


- (1 


-x)P+, 


Pci = x P 2 1 


P2± 


= (1 


- x)P 2± 





k'±,(20) 



which require that p^ = p^~ and pgj_ = p'q±. 

The hadronic matrix elements in the radiative decay 
between two particles A and X are of the form 



(X(P 2 )\J^\A(P l )) 

f 1 f d 2 k± 
= Aj / dx I — ^<f) 2 (x, k'±)<j) 1 (x,k ± ) 



16tt 3 



eEMK Q 2 )K(Pi ■ 1) - Pi(tv ■ ?)],(21) 



where k'j_ = k^ xq_±, Aj is the overlap of the j th 
flavor portion of the flavor wavefunctions and contains 
all of the relevant charge factors and mixing amplitudes. 
Comparing the last line of Eq. (|21() with Eq. I jlSjl it can 
be seen that we have defined the form factor in terms of 
the charge and mixing amplitude independent quantity 

liK,Q 2 ). 
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The sum of the light-front spinors over the helicities in 
Eq. H21|) is obtained as 



* x rr r+ 

\ Pi \ Pi 



Trf 



Pi'pi 



(22) 



The explicit form of the trace is summarized in the Ap- 
pendix. For V{P\) — > S(P 2 )^*(q) decay, the (transverse) 
polarization vector ey of the vector meson is given by 
ey(±) = [0,0, ej_(±)] and the trace term with the plus 
component of the current is given by 



S+= -{mAq R -^l-x)k R ]- 



2k 1 



x(M + 2m , 



[(2x - l) 2 m 2 + k 2 ± - xk ± ■ q ± 



where we use ey(+). Therefore, we obtain the one-loop integral, Ii(rrij,Q 2 ), as follows 



h(m j ,Q' 2 ) 



dx 



d 2 k 



16tt 3 

R„L\ /„2 



(^ s (a;,k'_ L )0y(j:,k_ L ) 



1 



(M + 2m,) [[2X L) m > 



x(l - x)M 
xk± ■ qj_] -+ 



xrrij[l — 4(1 — x) 



1 — x, k i 



k R q L 
-k±) 



r 



(23) 



(24) 



where k R q L = kj_ • qx — i\k± x q_L \z and even though the 
cross term does not contribute to the integral, the dot 
product term does contribute. 

Then, the transition form factor Pys 7 * (Q 2 ) is given by 

Fvsr'iQ 2 ) =A n h(m n ,Q 2 )+ A s h{m s ,Q 2 ), (25) 

where A n and A s are the overlaps of the up-down and 
strange portions of the flavor wavefunctions respectively. 



For example, in the case where V = <f> and S = (uu — 
dd)/y/2, An = -(sm.S u - ( j > )(e u - e d )/2 and A s = 0. Also, 
for V = 4> and S — ss, A n = and A s = — (cos5 w -^)e s . 

The transition form factor Fsv-y* {Q 2 ) f° r S(Pi) — > 
^(-^2)7* (9) can be obtained from Eq. (|24l) by replacing 
q^ — > q^ and kj^ — > k'^ and the explicit form for the 
one-loop integral corresponding to Eq. Ij24(l is given by 



J 



, f 1 f d 2 k± 1 f fc^g 1 ' 

Ii(mj,Q ) = I dx J j^cj)v(xy±)(t>s(x,k ± ) ^ i _^ M ^ xm j [-(2x-l) 2 -4(l-x)-^-} 



[(2x- l)W + ki - xk_L • qj + (x -> 1 - x; k -> -kj.) k (26) 



where Mq 2 = (m 2 + k']J/x(l - x) and again the cross 
term in k R q L does not contribute to the integral. As 
one may expect, however, we found that 7{(m,-,Q 2 ) = 
-h{mj,Q 2 ) and thus F S vr(Q 2 ) = -FvsyiG 1 )• 



B. The process S — > 77 

We now apply this model calculation to the two pho- 
ton decays of scalar mesons. In this case, the coupling 
constant gs-yy for S — > 77 is defined by 



M 2 = (7(ft)|ei • J|S(Pi)) 

= e 2 W(ei ' £ 2)(-Pi • q) - (d ■ Pi)(e a • ?)],(27) 

where ei = e 7 (<?) and £2 = e 7 (P2). In terms of <?s 77 , the 
decay width for this process is given by 

T = \a 2 g 2 ^Ml . (28) 

Here again, instead of calculating the two real photon 
decays, we first calculate the matrix element for S — > 77*, 
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which is given by and take the limit Q 2 — > to compute the decay rate 

for the two real photon decays. Using the same quark 
■^2 = (l(F'2)\J fl \S(Pi)), momentum variables in q + = frame as Eq. $Zj$ with 

= e 2 Fg 77 .((5 2 )[e^(P 1 • q) - Pf(e 2 • q)], (29) the plus component of the current, we then obtain 



m~2 = y^y^ij y^ idx 



Ai , A2 , A 



T JO 



16tt 3 



<fe(xi,kj_) 



V\(Pq) j U\ 2 (P2) U X2 (p 2 ) + U Xl (pi) 

<1 7= 



1 Pi VP2 



qi - [pL + m |]/P2 - [p|i. + 

^^■/aK, Q 2 )) [eftPt • «) - Pf (e 2 • <?)], 



— + (x — > 1 - kj_ — > — kj_) 



7?, 



(30) 



where A,- now contains factors of e 2 due to the presence of two electromagnetic vertices. As in the case of Eq. (|22H . 
the sum of the light-front spinors over the helicities in the first term of Eq. (|3U|I is obtained as 

TM = V\{Pg) ^ u\ 2 (p 2 ) u X2 (p 2 ) ^ + m Ai (Pi) 



ax,a 2 ,a ypg yp 2 yp 2 ypi 

Tr[(^ - mj) ^2(^2 + mj)7+(^i + m^)] 



2^2p±(p+)V 

4m 3 ■ e 2 - P2 ■ £2) + ggXgg -e 2 -pi- e 2 ) + pj"(p 2 • e 2 - pi ■ e 2 )] 

2 v /2pt(p+)V 



(31) 



where we have used the fact that (±) = 0. Now, using e 2 (+) = (0, V2q R / P^ , ex), we finally obtain the one loop 
integral, 



h{m,j,Q 2 ) = -a/6 dx — -4r(f>s(x, kj_) 



16tt 3 

(x — > 1 — x, kj_ — > — kj 



,[(2x - l) 2 + 4(1 - x)(fc V/qj)] x(l - x) 

(l-a;)Vx(l-x) m 2 + (k ± -xq ± ) 2 



(32) 



and the transition form form factor is given by 

F Sj t(Q 2 ) = A n I 2 (m n ,Q 2 ) + A s I 2 (m s ,Q 2 ). (33) 

As an example of the coefficients A n .s, consider the case 
where S = /q(1370). Here, A n = c 3 [(e 2 + e 2 )/\/2] and 
A s = b%e 2 , where C3 and 63 are the glueball-quarkonia 
mixing amplitudes of Eq. I|12|) . 



V. NUMERICAL RESULTS 

A. Decays involving /o(1370), /o(1500), and /o(1710) 

The expressions for the one-loop integrals, I\(mj,Q 2 ) 
and I 2 {m,j,Q 2 ), are evaluated numerically and used in 
Eqs. (|2"5|) and to compute the Q 2 -dependent tran- 
sition form factors for 77*, 4>j*, and pj* decays of the 



scalar mesons. As an example, we give the results for the 
case of the heavy glueball (i.e. M nfL < M sS < M gg ). The 
transition form factors for the 77* process are shown in 
Fig. [3 and those for the 4>j* and pj* processes are col- 
lectively shown in Fig. [21 In the case of the two photon 
decay, the form factor should fall off like l/Q 2 due to 
an intermediate quark propagator which becomes highly 
off-shell at large Q 2 . Figure 0] shows the behavior of 
Q 2 x P/ o77 * (Q 2 ) for each scalar state. Each of the curves 
clearly shows a tendency to flatten out, demonstrating 
l/Q 2 dependence in the form factors. 



The decay constants for the real photon processes can 
be obtained from the form factors in the limit as Q 2 — > 
(i.e. g — F(Q 2 — 0)). In this limit, the values of the 
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FIG. 2: fo — > 77* transition form factors for /q(1370) [dash- 
dotted], / (1500) [dashed], and / (1710) [solid]. 



FIG. 4: Q 2 times the fo — + 77* transition form factors (Fig.|5Jl 
for /o(1370) [dash-dotted], /o(1500) [dashed], and /o(1710) 
[solid]. 
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FIG. 3: fo — > P7* transition form factors for /o(1370) [long- 
dashed], / (1500) [dash-dot-dotted], and /o(1710) [short- 
dashed]; fo — > 07* transition form factors for /o(1370) [solid], 
/o(1500) [dash-dotted], and /o(1710) [dotted]. Here we have 
used cVtf, = +7.8°. 



one-loop integrals are 

/iK,Q 2 = 0) 
h(m s ,Q 2 = 0) 
h(m n ,Q 2 = 0) 
h(m s ,Q 2 = 0) 



2.05 GeV' 1 
1.93 GeV' 1 
-0.672 GeV- 1 
-0.375 GeV- 1 . 



(34) 



The decay constants are obtained by substituting these 



values into Eqs. I|25|l and (|33|l . The decay widths are 
then calculated using Eqs. i|17fl and l|28|). The widths for 
the 77, 07, and p7 decays for all of the glueball mass 
scenarios are summarized in Tables |HJ and IIIII respec- 
tively. The uncertainties in these values result solely from 
the uncertainties in the mixing amplitudes in Eqs. 1|13|) 
and l|14p. We have not accounted for the uncertainties 
in the meson masses. The uncertainties in the masses of 
/o(1500) and /o(1710) (~0.4%) are very small compared 
to the uncertainties in the mixing amplitudes (6%-40%), 
and can be neglected. However, the uncertainty in the 
mass of /o(1370) is about 10%, and would, therefore, 
contribute significantly to the uncertainties in the decay 
widths. 

Experimental data for radiative decays of the isoscalars 
/o(1370), / (1500), and / (1710) are poor. As one ex- 
ample, in the recent past the PDG had reported partial 
widths of 3.8 ±1.5 keV and 5.4 ± 2.3 keV for the process 
/o(1370) 77 24]. The PDG currently attributes these 
two values to /o(600), but at the same time they state in 
a footnote that this data could equally well be assigned 
to /q(1370) 1]. If these data which are on the order of a 
few keV do belong to /o(1370), this would be encourag- 
ing given that our results listed in Table I are consistent 
with this order of magnitude. However, the ambiguity 
noted above makes any such comparison irrelevant, and 
a great deal more experimental investigation is necessary 
before any definitive conclusions can be reached about 
the validity of any of the glueball mixing schemes. 

In the absence of good experimental data with which to 
compare our results, we turn to other theoretical predic- 
tions concerning these decay processes. InRef. 0]; Close 
and Kirk give predictions for ratios of fo — > 77 widths 
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TABLE I: Decay widths for the process fo — > 77. The unit 
of the decay width is [keV] . The uncertainties result from the 
uncertainties in the mixing amplitudes in Eqs. 11311 and 1141 . 





Light 


Medium 


Heavy 


/o(1370) 


1.6 


q+o-s 




/o(1500) 


8.0 


4 i+i-o 

*-l-0.9 


0.65lg;2 


/o(1710) 


0.92 


1 O+0.2 
1 -°-0.2 





which depend only on charge factors and mixing ampli- 
tudes, and ignore all mass-dependent effects. For the ra- 
tios r(/ (1710) -» 77):r(/ (1500) - T7 ):r(/ (1370) - 
77) they obtain 

Light Glueball = 1 : 5.1 : 2.8 
Medium Glueball = 1 : 2.4 : 3.6 
Heavy Glueball = 1 : 0.1 : 3.7 . (35) 

Our analysis which includes all of the relevant mass de- 
pendent effects yields 

Light Glueball = 1 : 8.7 : 1.7 
Medium Glueball = 1 : 3.2 : 3.0 
Heavy Glueball = 1 : 0.2 : 1.9 . (36) 

Our results differ slightly from those of Close and Kirk, 
however the same overall qualitative pattern is preserved. 

In Ref. HU , Close, Donnachie, and Kalashnikova 
(CDK) compute the 4> and p radiative decay widths for 
/o(1370), / (1500), and / (1710) in the NR quark model. 
Assuming that <p = ss, they obtain the values listed in 
Tables IIVI and Comparing these values with Tables 
ITT1 and [TTT1 it is clear that the relativistic corrections in- 
troduced by our model reduce the overall magnitudes of 
the decay widths by about 50-70%. We note somewhat 
greater reduction for the process of / (1370) — > <jry due 
to our non-zero 5^.^. A reduction in the widths would 
be expected given that the relativistic motion of the con- 
stituents tends to spread out the meson's wavefunction, 
thereby decreasing its peak value. Despite the differences 
in the overall magnitudes, however, the relative strengths 
between the different decay processes are fairly well pre- 
served. Just as in CDK's analysis, we find that the largest 
branching ratio is likely to be that of /o(1500) — > prf. In 
our model, this branching ratio is about 1% for the light 
glueball case, 0.6% for the medium glueball case, and 
0.2% for the heavy glueball case. 



B. Decays involving a (980) and /o(980) 

If we assume ao(980) to be a conventional qq, then 
the flavor structure should be (uu — dd)/y/2. For the 
processes ao(980) — > 77 and — > a (980)7, the decay 



TABLE II: Decay widths for the process fo — > <f)l- The unit of 
the decay width is [keV]. These values are for S u .^ = +7.8°. 
Using <5^_0 = —7.8° does not significantly alter the results. 
The uncertainties result from the uncertainties in the mixing 
amplitudes in Eqs. 1131 and 114H . 





Light 


Medium 


Heavy 


/o(1370) 


0.98 


0.831°;^ 




/o(1500) 


7.5 




170i^o 


/o(1710) 


450 


400+Jo 


36iH 



TABLE III: Decay widths for the process fo — > py. The unit 
of the decay width is [keV] . The uncertainties result from the 
uncertainties in the mixing amplitudes in Eqs. 113H and 114> . 





Light 


Medium 


Heavy 


/o(1370) 


150 


390±?o 


530±i?° 


/o(1500) 


1100 


630t«° 


210li 3 °o 


/o(1710) 


24 


trr + 16 
JO -14 


4ioif 6 ° 



constants and associated widths are calculated to be 

<?<Mo7 = -0.14 GeV-\ r 0ao7 = 2.8 eV 

5ao77 = -0.16 GeV-\ T ao77 = 990el/. (37) 

Our result for the magnitude of g^ao-y is consistent with 
the theoretical calculations of Gokalp and Yilmaz (2(|, 
who obtain 0.11±0.03GeV -1 using light-cone QCD sum 
rules, and Titov et al. 0, who obtain -0.16 GeV' 1 
from phenomenological considerations. However, none 
of these calculated values for the widths are consistent 
with experimental data. The <fi radiative width of 2.8 
eV gives a branching ratio of BR(<j) — > ao7) = 6.7 x 
10~ 7 which is significantly smaller than the PDG average 
of 0.88±£;^ x 10" 4 ; and, the two-photon width of 990 
eV is roughly 3 times larger than the value reported by 
Amsler of 0.30 ± 0.10 keV 0. The flavor content of 
the isoscalar /o(980) is less clear. If we consider the two 
possible extremes, /o(980) = nn and /o(980) = ss, we 



TABLE IV: CDK's results for f cjyy. The unit of the decay 
width is [keV] . 





Light 


Medium 


Heavy 


/o(1370) 


8 


9 


32 


/o(1500) 


9 


60 


454 


/o(1710) 


800 


718 


78 
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TABLE V: CDK's results for f 
width is [kcV]. 



P7. The unit of the decay 





Light 


Medium 


Heavy 


/o(1370) 


443 


1121 


1540 


Jo (1500) 


2519 


1458 


476 


/o(1710) 


42 


94 


705 



obtain 




= -0.05 GeV~ 


-l 




= 0.37 eV 


= -0.26 GeV- 


-l 




= 2.7 keV 


= +0.64 GeV~ 


-l 


r</>/o7 = 


= 60eV 


= -0.04 GeF- 


-l 


-^/o77 " 


= 63eV. ( 



ss 



The PDG average for the two photon width is T/ o77 = 
0.39^0^3 keV. Since the ss result falls below this value 
and the nn result sits above it, it would be possible for 
some mixed nn-ss state to reproduce the data. Work- 
ing out the mixing required for this, we find that /o(980) 
would be about 6% nn and 94% ss. This alone would 
allow for /o(980) to be interpreted as a conventional 
qq. However, the <fi radiative widths we calculated lead 
to the branching ratios BR(cj) — > /07) = 8.7 x 10~ 8 
for the nn and BR(<j) — > /07) = 1-4 x 10~ 5 for the 
ss. Both of these values fall well below the PDG av- 
erage of 3.3lg;5 x 10~ 4 . Also, if we compute the ratio 
BR((j) — > /o7) /BR((j) — > a 7) for our model predictions, 
we get 0.13 for f — nn and 21 for f = ss, while the 
experimental ratio is around 4. This again hints at the 
possibility of a mixed nn-ss being able to reproduce the 
data. However, the mixing needed to reproduce this ratio 
requires that /o(980) be 87% nn and 13% ss. This is the 
exact opposite of the mixing needed to reproduce the two 
photon width above. Overall our results are clearly in- 
consistent with the current experimental data on ao (980) 
and /o(980). Therefore, we conclude that these states are 
most likely not qq. 



VI. SUMMARY AND OUTLOOK 

We have performed the first LFQM calculations in- 
volving scalar mesons. First, the 3 P light-front wave- 
function was constructed. It was shown that, in gen- 
eral, the covariant operator used to obtain the spin-orbit 
wavefunction depends explicitly on the relative momen- 
tum between the meson's constituents, and is, therefore, 
more complicated than the naive form that is commonly 
used. This wavefunction was used to compute radiative 
decays involving / (1370), / (1500), /o(1710), /o(980), 
and ao (980). In the case of the three heavy isoscalars, 
the effects of glueball-gg mixing were taken into account. 
Specifically, three different mixing schemes correspond- 
ing to a heavy, medium, and light glueball were used. 



The lack of good experimental data made it difficult to 
draw any conclusions about which of the three mixing 
scenarios, if any, could be the correct one. We note, how- 
ever, that we have improve d up on the earlier NR model 
predictions of Close et al. [23. Relativistic corrections 
introduced by the LFQM resulted in decay widths that 
were about 50-70% smaller than those obtained in the 
NR calculations. Yet, very little change was observed 
in the pattern of relative strengths which is apparently 
quite robust. For the calculations involving ao(980) and 
/o(980), we assumed these states to be qq. In contrast 
to the case of the heavy scalars, there does exist well- 
established data for these light scalars. While one or two 
of the properties we caculated, when taken in isolation, 
could be considered consistent with the data, it is clear 
that our results as a whole do not match the data. This 
lends further support to the current idea that /o(980) 
and do (980) are not qq states. 

In a future work, we intend to refine our analysis and 
perform our own glueball-gg mixing calculation using the 
LFQM. In this current work, the model parameters (to, 
(3) used for the scalar ( 3 Po) meson wavefunction were ob- 
tained from a spectrum calculation fit to S-wave ( So and 
3 S±) meson data. In order to improve the scalar wave- 
function, we will perform a separate spectrum analysis 
using a QCD-inspired model Hamiltonian similar to that 
of Ref. [lj], and fit the spectrum to P-wave ( 3 Po, 3 Pi, 
3 P2, and 1 Pi) meson data. In addition to refining the 
model parameters, this analysis will also give the masses 
of the bare nn and ss P-wave quarkonia. With these 
masses, we will be able to perform a glueball-gg mix- 
ing analysis involving the isoscalars /o(1370), /o(1500), 
and /o(1710). These mixing amplitudes, obtained using 
the LFQM, could then be compared with those of Lee- 
Wcingarten and Close-Kirk. 



APPENDIX A: SPINOR STRUCTURE FOR 

V(S) -> S(V)j* TRANSITION 

In this appendix, we show the explict form of the trace 
given by Eq. J22J). For the V — ► S7* transition, the 
following two trace calculations are necessary 

s vsi = Tr[(pq-rn)(fl2 + rn)j + (tfi+m) ft 

= 4m[p^ (e ■ pg - e ■ p 2 ) + pt (t ■ Pq ~ £ ■ Pi) 

+Pq r ( e -Pl -t-P2) 

+ e + (pi ■ Pi + P2 ■ Pq - Pi ■ Pq ~ TO 2 )], 

S SV2 = Tr[(^ - m)(^ 2 + to)7+(^i +to)] 

= 4[pf {p 2 ■ Pq - m 2 ) + p^ [pi ■ p q - m 2 ) 



+pj(m 2 - Pl -p 2 )], 



to get 



S 



-1 



4(l-x)M 



Q + e • (Pi ~ Pq) a+ 
D VS1 ~ J ' 



M + 2m 



VS2 



(Al) 



(A2) 
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where e = e(-Pi) and we used the transverse polarizations 
in the calculation of the form factor and decay width. 

On the other hand, for the S — > U7* transitions, we 
have 



where e' — (.'{P2) and again we used the transverse po- 
larizations in the calculation of the form factor and decay 
width. 



9+ 
D svi 



S 



SV2 



to get 



Tr[(^g — m) / {i>2 + m)j + ($i + m)} 
4m[p^(e' -Pq~e' ■ p 2 ) + p 2 (^ ■ Pq - e' ■ Pi) 

+e' + (p 1 -p 2 +Pl -Pq~P2-Pq- TU 2 )], 

Tr[(^ — m)(p , 2 + m)j + (p' 1 + m)] 
4[pf(p2 -Pq - m 2 ) + p 2 (pi -Pq - m 2 ) 



+p g '(m -pi -p 2 )}, 



4(1 - x)M' 



D svi 



e' ' (P2 ~ Pq) c+ 



Mn + 2m 



SV2 



(A3) 



(A4) 



ACKNOWLEDGMENTS 

MD would like to acknowledge the support from the 
SURA/Jlab fellowship. This work was supported in part 
by a grant from the U.S. Department of Energy (DE- 
FG02-96ER 40947) and the National Science Foundation 
(INT-9906384) and Kyungpook National University Re- 
search Fund, 2003. The North Carolina Supercomputing 
Center and the National Energy Research Scientific Com- 
puter Center are also acknowledged for the computing 
time. 



[1] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 
[2] J. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 

(1982); Phys. Rev. D 41, 2236 (1990). 
[3] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985); 

R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987). 
[4] F. E. Barnes, Phys. Lett. B 165, 434 (1985). 
[5] C. Amsler and F. E. Close, Phys. Rev. D 53, 295 (1996). 
[6] V. V. Anisovich and A. V. Sarantsev, Phys. Lett. B 382, 

429 (1996). 

[7] R. L. Jaffe, Phys. Rev. D 15, 267 (1977); Phys. Rev. D 
15, 281 (1977); Phys. Rev. D 17, 1444 (1978); R. L. Jaffe 
and K. Johnson, Phys. Lett. B 60, 201 (1976). 

[8] F.E. Close, N. Isgur, and S. Kumano, Nucl. Phys. B 389, 
513 (1993). 

[9] N.A. Tornqvist and Matts Roos, Phys. Rev. Lett. 76, 
1575 (1996);N.A.T6rnqvist, Z. Phys. C68, 647(1995). 
[10] UKQCD Collaboration, G.Bali et al, Phys. Lett. B 309, 
378 (1993). 

[11] D. Weingarten, Nucl. Phys. B(Proc.Suppl.)34, 29(1994). 

[12] F. E. Close and N. A. Fornqvist, hep-ph/0204205 

[13] S. Nussinv and F. N. Fruong, Phys. Rev. Lett. 63, 1349 

(1989); Phys. Rev. Lett. 63, 2003 (1989). 
[14] H.-M. Choi and C.-R. Ji, Phys. Rev. D 59, 074015 (1999) 
[15] C.-R. Ji and H.-M. Choi, Heavy Ion Phys. 4, 369 (1996); 



Few Body Syst. Suppl. 10, 131 (1999); H.-M. Choi and 

C.-R. Ji, Phys. Rev. D 59, 034001 (1999); Phys. Lett. B 

460, 461 (1999). 
[16] Wolfgang Jaus, Phys. Rev. D 53, 1349 (1996). 
[17] W. Lee and D. Weingarten, Phys. Rev. D 61, 014015 

(1999). 

[18] F. E. Close and A. Kirk, Eur. Phys. J. C21, 531 (2001). 
[19] W.Jaus, Phys. Rev. D 41, 3394 (1990); Phys. Rev. D 44, 
2851 (1991). 

[20] P. L. Chung, F. Coester and W. N. Polyzou, Phys. Lett. 

B 205, 545 (1988). 
[21] C.-R. Ji, P, L. Chung, and S. R. Cotanch, Phys. Rev. D 

45, 4214 (1992). 
[22] H-.M.Choi and C.-R.Ji, Nucl. Phys. A 618, 291 (1997) 
[23] G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 

(1980). 

[24] D. E. Groom, et al, Eur. Phys. J. C15, 1 (2000). 

[25] F.E. Close, A. Donnachie, and Yu. S. Kalashnikova, 

hep-ph/02 10293 
[26] A. Gokalp and O. Yilmaz, Phys. Lett. B 525, 273 (2002). 
[27] A. I. Fitov, F.-S. H. Lee, H. Foki, and O. Streltsova, 

Phys. Rev. C 60, 035205 (1999). 
[28] C. Amsler, Rev. Mod. Phys. 70, 1293 (1998). 



